General Hubbard model for strongly interacting fermions in an optical lattice 

and its phase detection 
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Based on consideration of the system symmetry and its Hilbert space, we show that strongly inter- 
acting fermions in an optical lattice or superlattice can be generically described by a lattice resonance 
Hamiltonian. The latter can be mapped to a general Hubbard model with particle assisted tunneling 
rates. We investigate the model under population imbalance and show the attractive and the repul- 
sive models have the same complexity in phase diagram under the particle-hole mapping. Using this 
mapping, we propose an experimental method to detect possible exotic superfluid/magnetic phases 
for this system. 
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Among the control techniques for ultracold atoms, op- 
tical lattice and Feshbach resonance play particularly im- 
portant roles. The optical lattice is used to control the 
interaction configuration while the Feshbach resonance 
is a tool to tune the interaction magnitude. The com- 
bination of these two powerful techniques naturally be- 
comes the next frontier, which has attracted significant 
recent interest [l|, 0, H, 0, Q. To understand this im- 
portant system, one needs to have a Hamiltonian to de- 
scribe strongly interacting atoms in an optical lattice. 
The starting Hamiltonian is unfortunately complicated 
as one has to take into account multi-band populations 
as well as direct neighboring couplings 0, 0, M ■ We have 
described a method in Q to derive an effective lattice 
Hamiltonian for this system from field theory of the two- 
channel model. 

In this paper, we report the following advance along 
this direction: firstly, based on consideration of the sys- 
tem symmetry and its Hilbert space, we show that a lat- 
tice resonance model turns out to be a generic Hamilto- 
nian for this system. The resulting Hamiltonian agrees 
with the one from our previous microscopic derivation [j| , 
but the method used here shows this Hamiltonian should 
have general applicability. As an example, we point out 
that for strongly interacting fermions in optical superlat- 
ticcs, the effective Hamiltonian is again described by this 
lattice resonance model when we introduce some dressed 
degrees of freedom. For certain configurations of the su- 
perlattice, the system naturally supports d-wave super- 
fluid. Secondly, we mathematically map the lattice reso- 
nance Hamiltonian to a general Hubbard model (GHM) 
with particle assisted tunneling rates. The particle as- 
sisted tunneling brings in some new feature, in particu- 
lar, it may favor a superfluid phase compared with the 
Hubbard model. Thirdly, we discuss the attractive Hub- 
bard model with population imbalance between the two 
spin components, and show it has the same complexity 
in phase diagram as the repulsive Hubbard model under 
a particle-hole mapping. This result is related to the re- 
cent large effort to understand the polarized fermi gas 



Finally, using the mapping above, we propose an 
experimental scheme to detect possible exotic superfluid 
or magnetic orders in this system. The method is based 
on Raman-pulse-assisted time-of-flight imaging, and can 
reveal the superfluid or magnetic phases with detailed in- 
formation about the order parameter or the pairing wave 
function. 

For strongly interacting two-component (effectively 
spin-1/2) fermions in an optical lattice, when two atoms 
with different spins come to the same site, they form a 
dressed molecule with atomic population distributed over 
many lattice bands due to the strong on-site interaction 
[J, [5| . We consider the system with an average atom fill- 
ing number n < 2. In this case, we can neglect the 3- atom 
occupation of a single site as that is suppressed at low 
temperature by an energy cost about the lattice band gap 
Q. We then have only four possible configurations for 
each site i, either empty, or a spin-a (<r =|, J.) atom, or a 
dressed molecule. The creation operators for these con- 
figurations are denoted by b\ , a i(T , d\ , respectively, while 
the corresponding states are written as \b) i , ||, |) 4 , \d) i . 
We introduce the slave boson operator b\ for an empty 
site i so that the constraint of the Hilbert space on each 
site can be simply implemented through 
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Note that with this constraint, at a describe fermions 
while di and bi both represent hard-core bosons. 

We assume the system has a global SU(2) symmetry 
for the spin components. In that case, and |J,). 
are degenerate in energy, and the most general form 
of the single-site Hamiltonian can be written as Hi — 
a \a a i<? + (A- ~ where we have absorbed 

the single-atom energy into the definition of the chemi- 
cal potential /i, and A is the relative energy shift of the 
dressed molecule. For two neighboring sites i and j, due 
to the atomic tunneling and off-site interactions, there 
will be a Hamiltonian term to describe all the possi- 
ble configuration tunneling or couplings. With the spin 
SU(2) symmetry and the number conservation of each 
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spin component, the most general two-site Hamiltonian 
can be written as iJy = + where de- 

scribes the configuration tunneling that involves transfer 
of one atom with the following form (see the illustration 
in Fig. f ) 

H\f = ^(tal a b i b}a ja +t da dla ia at tT d :i J 

+ g(4 b 3 + djMKTOji - aa a it) + H - c -i ( 2 ) 

(2) 

and Hh ' describes the configuration coupling that in- 
volves real or virtual tunneling of two atoms with the 
general expression 

H if = \tdd\bib\dj + H.c.^j + x d n di n dj 

+ XaTiiHj + x s s l ■ Sj + x b n bi n hj . (3) 

(2) 

In , the number and the spin operators are de- 
fined by n dl = d\di, rii ee a^a^ + a^a^, n bi = b\bi, 

and s l ee J2crcr> <4 CT <W<W'/ 2 (<W is trie Pauli ma- 
trix). The term n^n^j is equivalent to the cross cou- 
pling n d irij + riin d j under the constraint (I). By analyz- 
ing the level configurations in Fig. 1, one can convince 
oneself that and include all the possible two- 
site coupling terms with the SU(2) symmetry. As the 
atomic interactions are short-range, all the multiple site 
couplings can be neglected. So a generic lattice Hamil- 
tonian is given by H = £. H t + (h$> + H^f) , 

where denotes neighboring sites. This Hamiltonian 
describes the coupling between the fermionic atoms ai a 
and the bosonic dressed molecules di with a detuning A, 
and will be referred in the following as the lattice reso- 
nance model. 




FIG. 1: Illustration of the configuration tunnelling between 
two neighing sites. The process shown in the figure correspond 
to the t, g, and td a terms in the Hamiltonian. 

The Hamiltonian H, together with the constraint (I), 
poses a well-defined problem. Note that H agrees in form 
with the effective lattice Hamiltonian for strongly inter- 
acting fermions that we derived before from a completely 
different method The only specification from that 
microscopic derivation is to fix the coefficients x b = 
and x a = —x s /4. As mentioned in in the case of a 
large detuning A, the Hamilton H is reduced to either 



the t-J model for atoms or the XXZ model for dressed 
molecules, depending on which species get populated. 
We also notice that for short range interactions, with 
increase of the lattice potential barrier, all the interac- 

(2) 

tion coefficients in decay much faster compared with 
those in ifj ■ . So for a lattice with sufficient depth, H$ 

(2) 

dominates over H-j, and in the following, without spe- 
cial mention we will consider the simplified Hamiltonian 

H = Ei Hi + E( W > b y dropping ffj». 

We now recast the Hamilton H into a different form 
which shows its connection with the Hubbard model. For 
this purpose, we map the dressed molecule state d\ \vac) 
to the two-fermion state a\^a\^ \vac), where \vac) denotes 
the vacuum. Note that physically the structure of the 
dressed molecules should be determined by diagonaliz- 
ing the on-site interaction Hamiltonian, and it generally 
involves superposition of atoms in many band configu- 
rations 0, [1| j which is certainly different from the state 
a iT a il \ vac ) with double occupation on a single band. But 
mathematically we can identify these two states by a one- 
to-one mapping. After this mapping, the Hamiltonian H 
can be written in the form 

H = [( A / 2 ) n * ( n * - !) - M n *l ( 4 ) 

i 

+ ^2 [t + Sg (n^ + rija) + Stn^n^] a\ a a 3<y + H.c. 

where Sg ee g — t, St = t da + t — 2g and riis = o,t~ai„ 
((j =1,1 for a —1,1). To verify the two forms of H in 
Eqs. (2) and (4) are equivalent to each other, one can 
check the physical process represented by each term to 
confirm it is identical. Note that in this new form of H, 
there is no need of the slave boson operator to constraint 
the Hilbert space as the latter is automatically fixed by 
the properties of fermions. As there is no additional con- 
straint, the Hamiltonian in the form of Eq. (4) looks 
simpler and may be easier for treatment in certain cases. 
Compared with the conventional Hubbard model, the ef- 
fective tunneling rate in H becomes an operator which 
depends on occupation of the two sites. The original lat- 
tice resonance Hamiltonian in Eq. (2) is thus mapped 
to a general Hubbard model ((GHM) with particle as- 
sisted tunneling rates. For weakly interacting fermions, 
the multi-band population and the direct neighboring 
coupling become negligible, then the coefficients g and 
t da tend to t a , and the GHM returns to the conventional 
Hubbard model as one expects in this case 0]. 

The derivation of the Hamiltonian H in this work is 
based on very general arguments about the single-site 
Hilbert space and the system symmetry. This reminds 
us that H has a generic form which should apply to dif- 
ferent systems with similar Hilbert space structure and 
symmetry properties. As an example, we point out that 
for interacting fermions in an optical superlattice, un- 
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der several interesting configurations, the system is also 
well described by the above Hamiltonian H . Figure 2A 
illustrates an optical superlattice potential which can 
be realized with two standing wave laser beams [Io| . 
With a combination of this superlattice and the con- 
ventional optical lattice potentials, one can realize the 
dimer or plaquette lattices as illustrated in Fig. 2B and 
2C where the intra-dimer (intra-plaquette) couplings are 
much stronger than the inter-dimer (inter-plaquette) cou- 
plings. To derive an effective Hamiltonian for this sys- 
tem, one needs to first construct dress energy levels for 
each dimer (plaquette) by exactly solving a few-site prob- 
lem. For two-component interacting fermions in those 
lattices near half filling, the low energy level configura- 
tions from each dimer (plaquette) have basically the same 
structure as those shown in Fig. 1 11, 12, IH, and the 
system also has the SU(2) symmetry. We then immedi- 
ately conclude that the Hamiltonian in the forms of Eq. 
(2) or (4) should be applicable to describe physics in the 
dimer or plaquette lattices around half filling. The two- 
dimensional plaquette lattice is particularly interesting: 
because of the internal plaquette structure, the excitation 
from \b) i to \d) i states in Eq. (2) has a <i-wave symmetry 

(e.g., \d\bij flips sign under a 7r/2 rotation of the lat- 
tice) [ll|, LL2J . When the effective detuning A is tuned to 
be negative, one can show under pretty well controlled 
approximations that the plaquette lattice in this config- 
uration supports d-wave superfluid [ll|, [13] ■ 



c -i 
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FIG. 2: Illustration of an optical superlattice: (A) The super- 
lattice potential. (B,C) The dimer and the plaquette lattices 
(bold lines represent stronger coupling) formed with the po- 
tential in (A). 

We now investigate some properties of the GHM in 
Eq. (4). When the detuning A is negative, similar to the 
Hubbard model, we expect this Hamiltonian is in a super- 
fluid state away from the unit filling. When A is positive, 
although it is not clear yet whether H has a superfluid 
state, compared with the corresponding repulsive Hub- 
bard model, we do expect that the superfluid possibility 
becomes higher when g > t a (which is very likely the 
case for fermionic atoms near a wide Feshbach resonance 
Q) as a large g term (see Eq. (2)) clearly favors Cooper 
pairing. From the single-site physics, we know that two 
atoms always have an on-site bound state (corresponding 
to a negative A) with the binding energy (—A) approach- 
ing zero as one moves to the BCS side of the Feshbach 



resonance 0, [Bj]. So in the ground-state configuration, 
the strongly interacting fermi gas naturally implement 
the GHM with a negative A. To experimentally investi- 
gate the GHM with a positive A, one needs to start with 
the population in atoms (instead of Feshbach molecules), 
and to approach the Feshbach resonance from the BEC 
side (the system is in a metastable state in this case). 
The effective Hamiltonians in different regions are shown 
in Fig. 3. 




FIG. 3: The effective Hamiltonians in different regions. The 
solid curved correspond to two dressed molecule bands, and 
the middle dashed line is an atomic band. On the BCS or 
BEC (with population mainly in atoms) sides, the Hamitoni- 
ans are given by the attractive (repulsive) general Hubbard 
models (AGHM and RGHM), respectively. As one increases 
the detuning A, one gets either the XXZ model for the dressed 
molecules or the t-J model for the atoms. In the deep BCS 
or BEC limit, the GHM returns to the attractive or repul- 
sive Hubbard model, and the XXZ Hamiltonian yields to the 
bosonic Hubbard model for molecules when multiple occupa- 
tion of a single site is allowed by weaker effective interaction. 

When we take into account possible population im- 
balance between the two spin components, the repulsive 
and the attractive GHMs (with positive or negative A, 
respectively) become intrinsically connected, and they 
should have the same complexity in phase diagram. Po- 
larized fermi gas recently raised a lot of interest [6( , and 
in free space (or in a weak trap) , although population im- 
balance yields some new features, the basic physics there 
is still largely captured by an extension of BCS type of 
mean-field theory Q. However, for polarized fermi gas 
in an optical lattice, we show that simple extensions of 
the BCS theory are very likely to give misleading results 
because of the exact mapping between the repulsive and 
the negative GHMs. Population imbalance corresponds 
to introduction of an effective magnetic field h, which 
adds a term —h'^2 li af (of = — mi) to the Hamilto- 
nian H in Eq. (4). We apply a particle-hole transforma- 
tion ajf — > ajf and an — > (— l)Vjj to the Hamiltonian 
[Tij ] (for simplicity, we consider a bi-partite lattice) . Un- 
der this transformation, the Hamiltonian H — h JT af is 
mapped to 



H' = J2h (A/2) m (m - 1) - n'af - tii 



(5) 



} J [ta- + 5g<j {nja + rija) + Stn^n^] a'^aja + H.c. 
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where the parameters fi' = fx — A/2, h! = h — A/ 2, 
t-f = t, t| = t^a) <5.9cr = g — ta, and we neglect the con- 
stant energy per site h — /j,. One can see that an attrac- 
tive GHM (A < 0) is mapped exactly to a dual repulsive 
model (with — A), where the chemical potential fi and the 
filed h exchange their roles. Superfluid phases (including 
both the BCS state and the LOFF (Larkin-Ovchinnikov- 
Fulde-Ferrel) state with pairing at nonzero momenta [151 ]) 
of the original model correspond to magnetic phases of 
the dual model and vice versa. For a repulsive Hub- 
bard model on a square lattice, the magnetic order exists 
only in a region near half filling, and with hole doping 
there are possibilities of exotic phases including a non- 
BCS superfluid state. This suggests for the attractive 
Hubbard model with population imbalance, the super- 
fluid phase exists in the region with small polarization. 
With further increase of the polarization, there could ap- 
pear exotic phases including a d-wave magnetic order. 
Experimental investigation of the attractive GHM with 
population imbalance (which might be easier for realiza- 
tion compared with the repulsive one) is therefore able to 
provide critical information to understand the challeng- 
ing phase diagram of the repulsive Hubbard model. 

Finally, we propose a method to detect possible exotic 
phases in this system by making use of the above map- 
ping. Our purpose is to directly measure the magnetic or 
superfluid order parameters. The detection scheme com- 
bines the time-of-flight imaging with some instantaneous 
Raman pulses [16(. We take 6 Li atoms as a typical ex- 
ample. The scheme is illustrated in Fig. 4. Right after 
turn-off of the trap, we immediately apply two consecu- 
tive impulsive Raman pulses. These pulses are assumed 
faster than the system dynamics (characterized by the 
Fermi energy), but slower compared with the level split- 
ting between the ||) and \[) levels (about 70 MHz). The 
first is a 7r-pulse, consisting of two laser beams propagat- 
ing along different directions, which transfers the atoms 
from the level ||) to |6) by imprinting a photon recoil 
momentum — q. As the level |6) is detuned from | j) by 
a few GHz, this transition at the same time tune the 
system out of Feshbach resonance (the atoms in states 
1 6) and |I) are only weakly interacting) The second 
is a 7r/2 pulse from co-propagating laser beams applied 
to the levels |6) and ||), which induces a transformation 
ak6 — * (a k 6 + e i¥; a ki )/\/2 and a k j. — * (a k | - e~ JV a k6 )/\/2 
that preserve the momentum k (ip is the relative laser 
phase). After these two pulses, we take the time-of-flight 
images (with basically ballistic expansion) for the atoms 
in levels |6) and ||), and the difference of these two im- 
ages give exactly the cross correlation of the f and J, 
spin-components at different momenta: 



?ik6 - ftkj. = 2 Re |e 1, '4+q,t fl kl 



(6) 



with this detection ability. (i) For magnetic phases 
with a pretty general form of the spin order parameter 
(sj) = V! cos(Q • r,j) + v 2 sin(Q • r,) we can con- 

firm it with sharp peaks for the correlation in Eq. (6) 
when the relative momentum q is scanned to ±Q . The 
spin vectors Vi and v 2 can be inferred from the rela- 
tive laser phase ip. (ii) For the LOFF superfluid state 
with pairing at a non-zero momentum q, the order pa- 
rameter (a k + qi |a_ k j) is nonzero. After the particle- hole 
mapping, this order parameter corresponds to a magnetic 
order 



,t 



-k-q,T l 



fl-kj, ) of the dual Hamiltonian. The peak 

of the correlation function in Eq. (6) at the relative mo- 
mentum — q thus confirms Bose condensation to a non- 
zero pair momentum for the original Hamiltonian, and 

the distribution in k of the correlation (|aL k _ q j»- k |^ 
gives the original pair wavefunction. (iii) Similar to the 
LOFF state, for a d-wave superfluid phase with the or- 
der parameter (a k) -[-a_ k j.) oc cosk x — cosfc y , , the pair 
wavefunction and its spatial symmetry can be directly 
measured by detecting the correlation (6) for the dual 
Hamiltonian. 



Pulse I 




Pulse II 




FIG. 4: Illustration of the two Raman pulses (a tt and a n/2 
pulses, respectively) before the time-of-flight for measuring 
the correlation function in Eq. (6). We use the level structure 
of 6 Li atoms as an example (with the magnetic field near the 
wide Feshbach resonance). 

In summary, we have established the results as we out- 
lined in the introduction. 

This work was supported by the MURI, the DARPA, 
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tracts, and the A. P. Sloan Fellowship. 



We now show through a few examples that we can 
directly confirm various magnetic or superfluid phases 
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